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Abstract 

We study vortex solutions in a holographic model of Herzog, Hartnoll, and 
Horowitz, with a vanishing external magnetic field on the boundary, as is appro- 
priate for vortices in a superfluid. We study the relevant length scales related 
to the vortices and how the charge density inside the vortex core behaves as a 
function of temperature or chemical potential. We extract a critical superfluid 
velocity from the vortex solutions, study how it behaves as a function of the 
temperature, and compare it to earlier studies and to the Landau criterion. 
We also comment on the possibility of a Berezinskii-Kosterlitz-Thouless vortex 
confinement-deconfinement transition. 
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1 Introduction 



One of the focal points of activity in applying holographic methods to condensed 
matter systems has been studies of holographic superconductors. An initial spark 
was the observation in [1], where it was noted that an Abelian Higgs model coupled 
to gravity in AdS space exhibits spontaneous symmetry breaking of the local U(l) 
symmetry (in the bulk) as the charge of the black hole was increased. In [2J, this 
phase transition was interpreted in the dual field theory as a superconducting phase 
transition, since in the model a charged operator gets a vacuum expectation value 



In this paper we take a more conservative stance and use the superfluid interpre- 
tation [Sill]. We are continuing our investigation of inhomogeneous extended config- 
urations in the model [2], which we begun in [51E] reporting on holographical dark 
solitons. In this paper we show that the model [2J allows vortex solutions even in the 
absence of a magnetic field on the boundary. In this case, the superfluid interpretation 
is most natural: in a superconductor a vortex would act as a source for a dynamical 
magnetic field and require turning on a flux line, whereas we are explicitly setting 
the boundary B = and yet find a consistent solution. Droplet and vortex solutions 



in the presence of boundary magnetic fields have previously been discussed in [?H9 



An early work on holographic vortices is a construction of a vortex line in pure AdS 
space [E] and in the AdS-Schwarzschild background [TB] . 

We construct vortex solutions in [2J, by solving the equations of motion of the 
Abelian Higgs model in the probe approximation, i.e. neglecting the backreaction 
of the scalar field and the gauge fields on the black hole geometry. We study basic 
properties of the vortex solutions, such as the free energy and associated length scales. 
We also study the amount of density depletion in the core of the vortex, and find 
agreement with earlier studies [5j|6]. In [6] we noted that the density depletion in 
the core of a dark soliton depends on the dimension of the condensing operator in 
a way reminiscent to a loosely bound fermion pair for dimension A = 2 and a more 
tightly bound fermion pair for A = 1, by comparing to an earlier study in non- 
relativistic superfluids [32]. For vortices a similar study of the depletion fraction in a 
non-relativistic setting is in |17j . 

In this work we extend the analogy between the type of the superfluid and the 
dimension of the condensing operator by studying critical superfluid velocities above 
which the superfluid flow starts to dissipate and superfluidity is ruined. The critical 
velocity is defined from the core size of the vortices. Our results for critical flows 
partially overlap with others obtained from constant superfluid flow in [3JII] . 

We also compare the critical velocity of superfluid flow to velocity of second sound 
(the second sound velocity was calculated in (3JHH] using different techniques), in the 

■'■For other studies of external magnetic fields on holographic superconductors see for example 



(VEV). 
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spirit of Landau criterion, and find that for both scaling dimensions the superfluid 
velocity is below the sound velocity. The proposal that superfluidity in the A = 2 
case might be due to loosely bound fermions is partly supported by the fact that the 
critical velocity is not set by the superfluid's second sound. 

The outline of the paper goes as follows. In section [2] we introduce the model 
and use cylindrical symmetry to derive the form of the vortex solution, and the 
corresponding equations of motion. In sections E] and H] we describe our numerical 
methods of solving the partial differential equations. In sections |5] and [6] we study 
the characteristic length scales of the vortex and the density depletion. Section [7] 
discusses the critical velocity as determined from the vortices. In section [8] we show 
that the free energy of one vortex is log divergent as is the case in superfluids (a 
related calculation has been done in [8]), and we study how the coefficient of the 
log term behaves as a function of temperature and comment on the possibility of a 
Berezinskii-Kosterlitz-Thouless (BKT) phase transition [19,20] in the system. Finally 
we end with a discussion and outline future directions. 

2 Review of superfluids and vortices 

We will begin by reviewing some basic properties of superfluids in 2+1 dimensions 
and how vortices arise in them and what role they play. A superfluid is characterized 
by spontaneous breaking of a global U(l) symmetry. This already poses a theoret- 
ical puzzle in 2 + 1 dimensions since such spontaneous symmetry breaking at finite 
temperature should not be possible due to the Mermin- Wagner theorem. In the holo- 
graphic model, this problem is perhaps being circumvented by having a large number 
iV of degrees of freedom in the field theory @. In this section, we will assume that 
spontaneous symmetry breaking is stabilized in the large- N limit and we will look at 
the consequences of that for the low energy physics of the superfluid. 

The low lying excitations in a superfluid are the Goldstone bosons from sponta- 
neous breaking of the U(l) symmetry. At large length scales, as compared to the 
inverse mass of the lightest massive quasiparticle, the system is well described by 
an effective action for the Goldstone field 4>(x). At zero chemical potential where 
the effective action is manifestly Lorentz invariant, the only non- irrelevant operator 
one can write down is (<9 M 0) 2 . Generalizing this to a finite chemical potential, fi, 
can be done by introducing a formal gauge field with one non- zero component, 
A = fj, [2S]. Since the effective action respects the formal gauge invariance — > e lA (f), 

2 See additional discussion in [2D - I23] . The question of identifying the N degrees of freedom can 
be addressed after embedding holographic model into string/M theory see for work in this 

direction. 
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— > + d^A it has the form 

S *ff = -J d 3 x^K(D^) 2 = j rf 3 x^((9 o 0-/i) 2 -(^0) 2 ), (1) 

The action (pQ) has a conserved momentum density T 0ai = Kdocpd^cp and a current 
density = Kdi<p. Both of the four vectors are parallel to d^, which motivates one 
to define the superfluid velocity vector v% = 7<9 M 0. The coefficients k and 7 seem at 
the moment to be arbitrary, but for a relativistic superfluid they can be connected to 
thermodynamic quantities as (30] 

«=^, 7 = -, (2) 
p p 

where p s is the superfluid density. These relations will be important in later sections, 
so that we will derive them in a simple example in the next section. 

We can define a vortex as a state for which the expectation value of the superfluid 
velocity (v l s ) = (d l cf)) / p is cylindrically symmetric and has a non-zero circulation 

c = p, J dxi(v % 8 ) , (3) 

where the integral is taken over circle of constant radius. Using the averaged (time 
independent) Heisenberg equation 

0=(d^) = fidi(v^, (4) 

we see that the radial part of the superfluid velocity is a constant. We will focus on 
the case where this constant vanishes. Using ([3]) we get the superfluid velocity of the 
vortex as 

(vl) = ( 5 ) 
P 

where p is the radial coordinate on two dimensions. For a normal fluid the circulation 
c can take any values, but in a superfluid (ft is the phase of a complex field and thus, 
the circulation is quantized in integer multiplets of 2tt. Quantized circulation is one 
of the hallmarks of superfluid vortices. 

Even though it is easy to see the existence of vortices simply by topological ar- 
guments and/or by use of the simple Goldstone effective action, the interior of the 
vortex is model dependent and carries interesting information about the microscopic 
structure of the superfluid [16J. Since the superfluid velocity diverges inside the vor- 
tex core, the vortex solution must interpolate all the way from the superfluid to the 
normal phase, with p s = 0, as a function of p. Thus, in order to describe a vortex, one 
is forced to appeal to a microscopic description that is able to describe the normal 
phase. 
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The energy of a single vortex can be estimated from the effective action ([T]) as 

E vor » J d 2 x\p Ps {{v s )f « ^log(y), (6) 

where _R C is an IR cutoff and £ a UV cutoff. In a microscopic description the UV 
divergence is absent since p s vanishes at the vortex core, so we are left with an IR 
divergence in the vortex energy. This does not still mean that vortices would not 
exist in superfluids. Vortices will be formed in a superffuid when the system is put 
under rotation with angular frequency u, which is larger than a critical frequency 
u c (see eg. [3 1J ) . If the size of the system is taken to infinity, the critical angular 
velocity tends to zero and vortices will immediately form when the system is rotating. 
Another scenario where vortices are important is during the BKT transition. In that 
situation, the entropy of a vortex becomes comparable to E vor /T and it becomes 
thermodynamically preferred to nucleate vortices in the system. 

The divergence of the vortex energy is in accordance with Derrick's theorem |32j . 
which states that the energy of a soliton solution in a theory with only scalar fields 
(with canonical kinetic terms) is divergent if the space dimension of the system D is 
larger than 1. In the holographic setting, there are also gauge fields involved, so that 
Derrick's theorem as such does not apply, and thus it is not a priori clear, whether 
the vortex energy is finite or not, in the holographic model. But as we will see, the 
vortex energy is indeed logarithmically divergent (for a related calculation see [8]). 
This is necessary for the superffuid phase to exist in non-zero temperature, since if it 
would take finite energy to create a vortex, it would be entropically favorable at any 
temperature to nucleate a soup of vortices, which would drive the order parameter to 
zero. 

2.1 An example model with vortices: The Gross-Pitaevskii 
equation 

To get some insight into the microscopic structure of superffuid vortices, before at- 
tacking the holographic problem, we study vortices in the relativistic Gross-Pitaevskii 
(GP) equation (we will drop the label "relativistic" in what follows). The GP equa- 
tion can be thought of as arising from a saddle point evaluation of a path integral for 
the bosonic field with the grand canonical action 

s GP = jd 3 x(\d ^- i p^\ 2 -m\ 2 -v\^\ 2 -^\ 4 y (7) 

The GP equation is relevant for weakly interacting, dilute Bose-Einstein (BEC) con- 
densates and tightly bound fermionic superfluids at low temperatures. 

It is interesting to see how one can end up on an effective action of the form 
(II]) from the GP theory. This gives a simple example where one can derive the 
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identifications ^ easily. By integrating out the fluctuations of the modulus in (J7J) 
we end up with an effective action for the Goldstone field , of the form 



s eff = J ^i^orfW-^i 2 -^, 



where \l/o is the VEV of This action is clearly of the form ([T]). In the GP equation 
the superfluid density is simply the particle number density derived from the action 
0, Ps = - ^d t ^* + 2ip J \^\ 2 ) w yu|^ | 2 , where in the second equality we 

have expanded around the ground state and ignored fluctuations. Thus the effective 
action becomes 

(9) 



S ef f = J rf 3 X^(|a o 0-/i| 2 -|^ 



|2 



which is the same thing as with k = p s j \i. The particle number current is now 
ji ~ |\l/o| 2, 9iX = Psdix/p = Ps v s, which tells us that v s = dtx/ P- Thus we see how the 
identification ([2]) arises from relativistic GP theory. 

Let us proceed to discuss vortices. For a time independent field configuration the 
GP equation is 

-V 2 ^ + (V - p, 2 )^ + g\^\ 2 ^ = 0. (10) 

Note that the non-relativistic version of the time independent GP equation can be 
obtained from fflUj) simply by replacing /i 2 by fi. An ansatz for a vortex solution can 
be taken in the form \1/ = e m0 R(p), where n is the winding number (or circulation) 
of the vortex. Plugging the ansatz into the GP equation gives 

1 n 2 

- R" R' + —R + (V — fi 2 )R + gR 3 = 0. (11) 

P P 2 

Again the superfluid velocity has the same form as in ([5]). Since ffTTl) does not have 
known analytic solutions for non-zero winding n, we will solve it numerically. The 
numerical solution is displayed in Fig [TJ From the numerics it seems that the field 
R(r) behaves in a power law fashion at large values of p. This suggests that we can 
try to get an analytic solution for asymptotically in p. Plugging a power series ansatz 
to (TTTj) gives the leading behavior 



and the particle number density (which is equivalent to the superfluid density) is 
given by 

,,2 l - ,,2, 

9 v P 2 



Ps { P ) = »R\p) = p^—^(l ~^$ + ...), (13) 



where £2 = 1/ y p 2 — V. For the unit winding vortex we can also define a length scale 
from the vortex's core, R ~ _R(oo)(p/£i + ...). We can numerically calculate £1, which 
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Figure 1: Spatial profile of a vortex solution to the GP equation. 



gives the ratio 



6 
6 



1.75. 



(14) 



Note that there is really only one length scale in the problem. This means that at 
sufficiently low temperatures (that is, when the GP equation is valid) the ratio ( |T4"|) 
is a constant as the chemical potential is varied. We will revisit this in the context 
of holographic superfluids in section [5j 



3 The model 

As argued in [2] , a holographic dual of a superffuid is provided by a scalar field coupled 
to an Einstein-Maxwell system in asymptotically AdS space. 



SAdS 



d A x\J—g 



-if UV*I 2 + m 2 l*| 2 + \f, 



L 2 



(15) 



M \ - / T 

As in [2J, will choose the mass m 2 = —2/L 2 . Throughout this work we will work 
in the probe limit, (-4 is small), so that the backreaction to gravity can be ignored. 
Thus, we are working with the fixed bulk metric 

L 2 



^r(-f(z)dt 2 + f{zY Y dz l + dp 1 + p 2 d6 2 ' 

z 2 



where f(z) = 1 — The AdS-CFT dictionary tells us that 

-S A ds(onshell) _ I - J d 3 x(p(x)fi(x)+Oi (x)Oi (x)) \ 

\ c IQFTi 



(16) 



(17) 



where the boundary quantum field theory operators are related to the boundary values 
of the bulk AdS fields by the relations 

V(x,z) = zOi(x) +z 2 2 (x) + ... A t (x, z) = n{x) + zp(x) + ... (18) 
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Thus, in order to find the QFT operator expectation values we need to solve the 
classical equations of motion in AdS space to obtain the on shell fields. Ignoring 
gravitational backreaction, the equations of motion become 

= -^d p {^g~g^d u R) + m 2 R - R(d pX ~ A p ) 2 (19) 



-9 

= -L=d^^Fn - R\A V - d v X ) (20) 

= d lt (y/=gieg^{d v x-A v )) i (21) 

where we have defined the real valued fields R and x according to the relation \I> = 
±Re% 

The equations of motion may be greatly simplified using the observation that 
cylindrical symmetry implies that gauge invariant quantities are independent of 9. A 
detailed discussion of gauge fixing and the use of symmetry is presented in Appendix 
lAl An important point discussed in Appendix |B] is that the field % is a normalizable 
mode. If x was non-normalizable, its boundary value would be set by hand, which 
in turn would imply that features such as the superfluid velocity profile are not 
determined by the system dynamically. 

Working in the A z = gauge and defining the field R = R/z, the equations which 
describe a vortex profile are 

= fd 2 R + d z fd z R -zR+ -d p (pd p R) 



P 

I 2 (Ae-nf 



~ R{~ L A + (22) 



= fd 2 z A t + -d p (pd p A t ) - R 2 A t (23) 
P 



= d z (fd z A e ) + pd p (-d p A e ) - R 2 (A e - n), (24) 

where all the fields are functions of only z and p. As discussed in appendix |A] the 
only non-zero gauge field components are A t and Ag. § 



3 Since there is a magnetic flux f F ^ going through the vortex core in the bulk, one can ask 
what happens to it as one approaches the boundary. Since the magnetic flux is conserved J s F = 
for any closed surface E, the magnetic flux has to go somewhere as one approaches the boundary. 
Since we have a vanishing boundary magnetic field the magnetic flux does not reach the boundary, 
but rather it turns to the p direction in the form of the magnetic field B p ~ In the holographic 
dictionary B p = d z Ag/p is not identified as a magnetic field, but rather with the superfluid current 
je. So, in short, the magnetic flux turns into superfluid flow on the boundary. 



8 



4 Finding the vortex solutions 



4.1 The method 

The equations of motion are a set of coupled nonlinear partial differential equations, 
which do not seem to be solvable analytically. In order to proceed, we therefore resort 
to numerical methods. Because the vortex solutions are inhomogeneous, it is difficult 
to use standard differential equation solvers, such as Mathematica's NDSolve. Instead 
we will use a Gauss-Seidel relaxation scheme. As described in 00, we first place 
the system on a finite box of radius L, (z,p) € [0,1] x [0, L] and discretize along 
both axes. After the usual procedure of replacing derivatives with finite difference 
derivatives, the differential equations are turned into finite difference equations. 

The resulting finite difference equations are solved by choosing an initial seed 
configuration and iterating with the Gauss-Seidel method. The error in solving the 
equations of motion falls off with a power law as the lattice size is increased. With a 
sufficient number of iterations and a fine enough lattice, the seed configuration relaxes 
to a solution to the equations of motion. For further details on this numerical method 
sec [6J. 

4.2 Boundary conditions 

To find a unique solution to the finite difference equations, we specify boundary 
conditions at each edge of the region being simulated. The boundary conditions at 
the AdS boundary determine the external sources turned on in the dual field theory. 
The gauge field, A t , takes a constant value (independent of p and 6) on the boundary 
corresponding to a constant chemical potential, A t (z = 0) = p. The scalar field, 
R = |^|/ '(\/2z), can satisfy either Dirichlet R(z = 0) = or Neumann d z R(z = 0) = 
boundary conditions, corresponding to two different boundary theories with different 
scaling dimensions of the condensing scalar operator, either (9 2 or 0\. Finally the 
angular gauge field, Aq, satisfies Dirichlet boundary conditions, Ag(z = 0) = 0. This 
is because we want the boundary magnetic field to vanish, implying d p Ag(z = 0) = 0. 
Because the boundary value of A e is simply a constant, we can choose such a gauge 
that A e (z = 0) = 0. 

On the p = L boundary of the region being simulated we impose Neumann bound- 
ary conditions for all the fields, since we want the fields to asymptote to the trans- 
lationally invariant symmetry breaking solution. We have verified that the solutions 
and results in this paper are unchanged when one increases L, hence we can safely 
conclude that we have used a large enough box. 

4 In Appendix [C] we discuss a complementary approach valid away from the vortex core. In that 
regime, the differential equations reduce to ordinary differential equations and one may use "off the 
shelf' tools. 
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Since the core of the vortex should be in the normal phase we impose Dirichlet 
boundary conditions on the scalar field R(p = 0) = in the core. Because the 
scalar field vanishes in the origin, there is no bulk electric charge at the vortex core. 
Therefore, along the p = boundary we impose Neumann boundary conditions for the 
gauge field A t such that center of the vortex is not a source for the radial electric field. 
This is basically a regularity condition. Finally we impose the Dirichlet boundary 
condition A e (p = 0) = on the angular gauge field. This may also be regarded 
as a regularity conditions along the vortex core. We want the bulk magnetic field 
B z = dpAg/p to be finite at the origin. This means that Ag has the behavior Ag = 
A ( 'q\z) + p 2 A^\z) near the vortex core. To have finite superfluid current in the core 
we require the d z Ag/p to be finite in the limit p — > 0. This means the d z A^ = 0. 
Furthermore, since we required Ag(z = 0) = fixes the constant part A d = and, 
thus A e (p = 0) = 0. 




Figure 2: Typical field profiles for the < 0\ > condensate. (For visualization purposes 
we have subtracted a linear background from A.) 




Figure 3: Typical field profiles for the < O2 > condensate. (For visualization purposes 
we have subtracted a linear background from A.) 

The final boundary conditions to discuss are regularity conditions to be imposed 
along the black hole horizon. The equations of motion are elliptic outside the horizon, 
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but become parabolic along the horizon. Also, we require that A t (z = 1) = so that 
the connection one form be normalizable at the horizon. This leads to the following 
regularity conditions at the horizon. 



o = (dj)d z R -R+ -d p ( P d p R) - R K -^- 

= A t 

= (d z f)d z A e + pd p (-d p A e ) - R 2 (A e - n). 




(25) 



(26) 



(27) 



4.3 Existence and stability 

Vortex solutions of (|22|) -( |24|) are topologically stable, with a conserved topological 
charge 



For all regular solutions n is an integer, and thus continuous time evolution cannot 
change its value. 

If one prepares the system in a configuration of non-zero winding, the time evolu- 
tion will keep the field in the same winding sector for all times. Eventually the system 
will stabilize to the minimum energy field configuration which solves the equations 
of motion with vanishing time derivatives. They are simply the static equations 
(l22l) . ( !23l) . and ( )24l) . The minimum energy field configuration is approximated by the 
solutions we find numerically. 

We are solving the equations using a relaxation method, which means that the 
solution evolves in the iteration time r by a (generalized) diffusion equation. In order 
to establish that our numerical solutions approach the minimum energy solution in 
the desired winding sector, we need to show that the field R does not generate a 
discontinuity to "unwind" the solution. The fact that this does not happen follows 
from the diffusive nature of the iteration algorithm. The diffusive nature guarantees 
that the energy of a configuration is non-increasing. Thus, if we start from a seed 
configuration with finite energy density above the ground state, we will not generate 
discontinuities to the fields since they would require too much energy. 

5 Characteristic scales 

Typical condensate, charge density, and superfluid density radial profiles as extracted 
from the bulk solutions, are shown in Fig.(J3|). 

We begin our analysis of the vortices described in section (14. ip by focusing on unit 
winding vortices. The first quantities we would like to identify in these vortices are 
the relevant scales of variation. It can be seen from our numerical solutions that the 
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fields have a power-law fall off as p — > oo. This fact is confirmed in the tail regime as 
is discussed in appendix O 

Using the standard AdS/CFT dictionary, the behavior of the condensate may be 
read off from the boundary behavior of the scalar field. In this way we can define a 
coherence length by either fitting a power-law function of the form 

(O i (p)) = <P\(O i {oo))\{l-& + ...), (29) 

to the tail of the condensate field or from the slope of the condensate in the vortex 
core. 

(O t (p))=e»\{O t (oo))\£ + ... (30) 

As discussed in Appendix [UJ in the large p regime the differential equations may 
also be solved using shooting methods. The length scale in ( 1291) may be determined 
with high accuracy. We will conventionally define the length scales as was done for 
solutions of the GP equation 12.11 

The boundary value of the derivative of the gauge field A t is identified with the 
boundary theory charge density p q (or the total particle number density in the super- 
fluid language). It is seen that it can be fitted well to a function of the form 

p q (p) = PM + Mp) - P,(0))(1 - f 2 + ...), (31) 
with a characteristic charge density coherence length £ 9 . 



<o> <P?> K P s 




Figure 4: Typical radial profiles for expectation values. 



The boundary value of the derivative of the angular gauge field Aq can be identified 
with a current ]q in the boundary theory. We will define the superfluid density p s 
through the relation 

f = Psvl (32) 
Since x( z — 0) is the phase of the operator expectation value 

(Oi) =Rie ix{z=0) , (33) 
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Figure 5: The length scales related to different operator expectation values. The 
different colors blue (up), black (middle) and red (down) correspond to the scales of 
the condensate, superfiuid density and charge density, respectively Left is for 0\ and 
right is for 02- 



we can identify the superfiuid velocity as 

vl = -d l x(z = 0) = -^^, (34) 
fi P P 

where the normalization factor follows from relativistic symmetry [30J as we men- 
tioned earlier. Now the superfiuid density for the vortex can be identified as 

pd z A e (z = 0) 

ps = • (35) 

n 

We can identify a third scale by fitting the superfiuid density as 

p.(p)=p a (oo)(l -§ + •••)■ (36) 
P 2 

In Fig. ([5]), the different scales as determined from the asymptotic forms fl29|) . fl3Tl) 
and ( l36i) are graphed. As was found for holographic dark solitons, the length scales 
typically have different dependences on the chemical potential, indicating that there 
are independent scales governing the physics. There is one exception to this result, 
for the 02 superfiuid, the condensate and charge density length scales coincide. 

Another interesting quantity is the ratio of the condensate length scale determined 
from the core to the one determined from the tail £i/£2- It was noted in [T7j that for 
a BCS superfiuid, near zero temperature, these two length scales are very different. 
The basic idea behind this is that the physics at the vortex core is determined by the 
normal phase properties of the superfiuid, so that the length scale is related to the 
microscopic Fermi momentum, £i ~ l/kp, while the length scale in the vortex tail 
is determined by the physics of the superfiuid phase, where the relevant scale is the 
inverse gap £ 2 ~ 1/A. Thus, in weakly coupled BCS theory £x <C £ 2 - F° r a more 
tightly bound superfiuid it was found in [T7j that the two length scales coincide at 
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Figure 6: The ratio of the core length scale to the tail length scale of the condensate 
as a function of the chemical potential. 0\ is blue and O2 is red. 

the unitarity limit (infinite scattering length) and on the BEC side of the BCS-BEC 
crossover Jf] 

Our results for the ratio ^1/^2 i n a holographic superfluid are shown in Figure 
|6j There are two features worth noting. First, for the 0\ superfluid the ratio of 
tail to core length scales is consistent with a constant ratio at low temperatures 
(large fi). As for the GP equation, we estimate the core length scale using the 
slope of the condensate at the origin. There is a relatively large uncertainty in the 
finding the slope in the core because the spatial gradients are large (especially at 
large fx). The uncertainty in the length scales may be estimated as in [6]. We find 
that the uncertainty ranges from 1 — 10% with increasing chemical potential. For 
the O2 superfluid the ratio of length scales clearly depends on the chemical potential, 
consistent with there being two distinct length scales in the vortex's core and tail. 

The second feature we note is that the ratio of the length scales is larger for 
O2 than for 0\. Given that other quantities (including the density deplitions and 
the critical velocities) indicated that O2 might be a BCS-like superfluid. This is a 
surprising feature. In a speculative vein, if we assume that the 2 liquid is a fermionic 
superfluid, we would interpret the absence of a small core length scale as suggesting 
a vanishing Fermi momentum. This is also supported by the fact that we do not see 
any Friedel oscillations in our solutions, which would have a wavelength proportional 
to hp |lb[ITT]. One possibility is that, there are a large number of charged fermion 
species N contributing to the condensate. In this way, one could put an infinite 
charge to the ground state without having to occupy higher energy states (assuming 
the charge density is finite), which would mean that the Fermi surface, and thus the 
Fermi momentum would scale to zero as N — > 00. 

5 We found in section 12.11 that the ratio of the two length scales is not one, but this is not really 
to be expected, but rather both of them are proportional to the same length scale £ = 1/vV 2 — V. 
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Figure 7: The density depletion fraction as a function of temperature. The upper 
curve is 0\ (blue) while the lower curve is O2 (red). 

6 Density depletion 

An interesting observable in the vortex solutions is the charge density profile. It carries 
non-trivial information about the fraction of the charged matter in the condensate. 
Near T = T c the density depletion in the core of the vortex is very small, which 
suggests that the part of charged matter in the condensate is small compared to the 
total charge density. This is found for both condensing operators 0\ and O2. 

A pronounced difference between 0\ and O2 becomes clear when the temperature 
is lowered. For O2 the density depletion fraction seems to be saturating at 40% while 
for 0\ it is likely to grow near 100%. This same basic pattern was also observed for 
holographic dark solitons |5]|6|. Comparing the density depletions to those obtained 
in a non-relativistic setting [T7] suggests that one may identify 2 as a BCS type 
superfluid and 0\ as a BEC type superfluid. 

7 Critical velocity 

An interesting quantity in the context of superfluids is the critical superfluid velocity, 
above which the superfluidity of the system is destroyed and the flow of the fluid 
starts to dissipate. 

7.1 Landau criterion 

At low temperature the critical superfluid velocity can be estimated by the Landau 
criterion, which goes as follows. Consider the superfluid moving in a container with 
a velocity v. Dissipation occurs when the wall of the container or some other defect 
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Figure 8: The left figure shows a typical profile for the supercurrent in a vortex. The 
right figure shows the critical radius for different values of n and it seems to be a 
linear function. 



starts to excite quasiparticles in the fluiqj. Let us consider a quasiparticle with 
energy e p and spatial momentum p. Obviously these together form the momentum 
four vector. As we boost back to the rest frame of the container, the time component 
of the four momentum transforms into 



2 



(37) 



Whenever the energy of the quasiparticle e' is negative, as measured from the rest 
frame of the container, it becomes energetically favorable to create such excitations. 
The expression ( |371) is minimized when p and v are antiparallel. This gives us the 
Landau critical velocity 

v crit = nnnp-rjp (38) 

where the minimum is taken over all possible quasiparticle excitations. This shows 
that it is not necessary for the quasiparticle spectrum to be gapped, but there can be 
excitations with linear dispersion relation for small momenta. 



7.2 Determination of the critical velocity from the vortices 

We can determine a critical velocity from the vortices as follows. We begin by asking 
a simple question: Why does the radius of a vortex increase as the winding number 
is increased? The simplest physical reason is that the velocity of the superfluid flow 
around the vortex increases as v ~ n/p and at some radius p* the superfluid velocity 
gets larger than the critical velocity and inside that radius the condensate vanishes 
since it is no longer energetically favorable. 

Within our numerical solutions we can easily test this idea, by plotting the critical 
radius p*(n) and seeing whether it behaves as a linear function of the winding n, as 
should be if the idea of critical velocity is the correct physical reason for the grow of the 

6 For the dissipation one does not need a container, it is enough to have relative motion between 
superfluid and normal components. 
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Figure 9: The upper curve is the sound velocity and the lower curve is the critical 
velocity as a function of the temperature as determined from the vortices. 0\ is left 
and O2 is right. 

vortices. Following ideas from |17] . we define the critical radius to be the point where 
the current j l (p) reaches its maximum. A typical profile of the supercurrent j s (p) for a 
vortex is shown in Fig.®. Inside the a critical radius, where the supercurrent reaches 
its maximum value, the superfluid density starts to decrease causing the supercurrent 
to decrease. The critical radius as a function of the winding number, p*(n) seems to 
be indeed a linear function of n as shown in Fig.®. 

By generating vortex solutions with a high winding number n = 20, for different 
values of the chemical potential (or equivalently, for different values of the temper- 
ature) we can see the behavior of the critical velocity as a function of the chemical 
potential. The results are shown in Fig.(|9]) for different condensing operators 0\ and 



Near the critical temperature the functional form of the critical velocity fits well 
with a square root behavior 



This is the same functional form as was found in [3J . It is interesting to compare the 
critical velocity determined by the vortices to the Landau critical velocity. For bosonic 
superfluids and tightly bound fermionic superfluids, the quasiparticle determining the 
Landau critical velocity is a collective mode (sound mode) [33]. A vortex's critical 
velocity depends on the type of superfluid. For a weakly coupled fermionic superfluid 
it is the breaking of Cooper pairs that determines the critical velocity, while for a 
bosonic superfluid a sound mode sets the relevant scale. In [31] it is found that 
the second sound is the smallest (known) sound velocity. Thus, it should determine 
the Landau critical velocity for a bosonic superfluid. We can use this to probe the 
superfluid type. In Fig. ([9]), the critical velocities are plotted with the corresponding 
sound velocities. 

At low temperatures (large /z), for the 0\ superfluid, the critical velocity is closer 



o 2 . 




(39) 
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to the sound velocity than for the O2 superfluid. This indicates that the vortex critical 
velocity is more likely to be set by a sound mode for the 0\ than for (9 2 - 



8 Free energy and a puzzle with the BKT argu- 
ment 

Next, we will show that the vortex solutions have a logarithmically divergent free 
energy as is usual in superfluids (for a related calculation see jH]). According to the 
AdS/CFT dictionary the free energy of the boundary QFT is identified with the 
Euclidean on shell action in the bulk AdS. First we will evaluate the Lorenzian on- 
shell action on the vortex solution, with a cutoff z = e at the AdS boundary and an 
infrared cutoff at the radial position p = R c . 

s AdS = J dt J\z J\e £ c P dp(- ^\d z v\ 2 - -^\d p ^\ 2 



+ — {d p A t f - ^{d z A e f - -^\d e V - iA e ^\ 2 ) (40) 



After substracting a counterterm S ct = J d 3 x\^f\ 2 /e 3 the Lagrangian density is finite 
everywhere for the vortex solutions, but the action diverges due to integral over p. 
Thus, leading terms in the free energy are the diverging ones as R c — > 00. We will 
concentrate here on these terms. Because the diverging terms are a large p effect, 
they may be captured using the asymptotic expansion in Appendix O 

* = + A« = 4»W + ^#, A, = A°(z) + S -M?l. (41) 

pZ pA pZ 

At leading order in this expansion the Lagrangian density is simply £(R = R°,A t = 
A®,Ae = 0). This gives rise to an extensive term in the free energy diverging as R 2 , 
which is exactly the free energy of the translationally invariant symmetry breaking 
state. At subleading order there are logarithmic divergences in the free energy. These 
have two sources, the two last terms in (l4"0l) are log divergent when evaluated on 
Ag(z), while the power law correction to Ag in (14 ip gives finite subleading terms. 
Another source of possible log divergences are the first order corrections to R and A t 
in (JH]). It is easy to see that these terms have to vanish, since they give rise to terms 
that are proportional to the asymptotic equations of motion (1731) . schematically 

\5A t p* + 5 R P 2 r ^ 
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Figure 10: Free-energy coefficients as functions of the chemical potential. Blue (up) 
is 0\ and red (down) is C>2- 



and to boundary terms 5A(0)d z A^(0) and 8R(0)d z R(0), which both are vanishing. 
Thus, we see that the difference between action evaluated on the translationally in- 
variant symmetry breaking solution and the vortex is given by 

AS = -27rlog(^) J dtdz^-{R°)\n - A° e ) 2 + f -{d z A" e f] , (43) 

where we have neglected terms that are finite in the limit p — > oo. By using the 
equations of motion of the Ag field we see that the integrand in (|4"3"1) becomes simply 
nd z (fd z Ag). By continuing the action to Euclidean time with period /3 = 1/T we end 
up with the free energy difference 

Afi = «log(^), (44) 

where a = irn(d z Ag)\ z=0 . The coefficient a is shown in figure [TU] as a function of the 
chemical potential. In the section [5] we showed that the superffuid density for the 
vortex solution is simply p s = p(d z Ao)/n. Thus, we see that the vortex free energy 
is given by 



An 



Tip s (oc)n 2 /R, 



log(f), (45) 



where £ is a length scale which measures the vortex core size. The vortex free energy 
agrees exactly with the one determined from symmetry arguments in section [2j as it 
should. Furthermore, a/27r is the coupling constant in front of the Goldstone effective 
action for the holographic model. 

There is a slight puzzle involved with the free energy as it vanishes at T = T c . Even 
though we have calculated the vortex free energy, we are missing a piece involving 
the vortex entropy, since we have fixed the center of the vortex. The entropy of a 
single vortex is simply the number of possible vortex states. Since the vortex occupies 
an area of order £ 2 , the number of possible vortex positions is R 2 C /^ 2 and thus, the 
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entropy is S = log(-R^/£ 2 )[] The free energy of a single vortex is 

F = E -TS = (a-2T)log(y). (46) 

This suggests that at some temperature below the critical temperature of the phase 
transition described in [I], it becomes entropically favorable for a Berezinskii-Kosterlitz- 
Thouless (BKT) transition. In such transitions it is entropically favorable for vortex/ anti- 
vortex pairs to deconfine, destroying the superfluid state. A potential pitfall of the 
above argument is that the semi-classical gravity approximation involves a large- N 
limit. There is the possibility of an order of limits issue when taking both N and 
q — > oo. It is possible that factors of N enter into the ratio E/T but do not effect 
the entropy Su In this case there would be an order of limits issue, which might 
drive the BKT transition temperature to the superfluid transition temperature. To 
clarify this issue it may be necessary to include 1/N corrections. We will save a more 
detailed study for later work. It would also be interesting to see where the potential 
BKT transition is placed in the phase diagrams of [35] . 

9 Discussion 

In this work we have studied vortex configurations permitted in holographic super- 
fluids. These vortex solutions are characterized by a depletion of the condensate 
and charge density in the vortex core and the quantization of the circulation around 
the vortex. In the bulk these solutions correspond to vortex configurations in an 
Einstein-Maxwell-Higgs system, which have vanishing boundary values for the spa- 
tial components of the vector field. Bulk fluxes do not become boundary gauge fields, 
but instead map into the charge and superfluid densities in the boundary theory. 

We explored features of the vortices as functions of the chemical potential and 
type of condensing operator. We have used the fact that vortex configurations are 
not small deviations from the homogeneous condensate to probe both long and short 
distance features of holographic superfluids. There has been much work interpreting 
the system of [2] as a superconductor. This may be appropriate for features which 
do not require a dynamical gauge field (such as conductivity calculations and the 
superconducting gap). However, for other features, such as the Meissner effect, a 
dynamical gauge field and Maxwell's equations are required. Then, vortices act as 
sources for the gauge field and require magnetic fluxes. However, the existence of 
vortex solutions even in the absence of a magnetic field on the boundary indicates 

7 Rather than appeal to this heuristic derivation of the entropy one can see the same term arising 
when one path integrates the moduli associated with translating the vortex core in the p — 6 plane. 

8 S is determined by the metric on moduli space and its leading contribution is independent of 
coupling constants in the gravitational theory 
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that these systems may more likely be a superfluid or a thin film superconducting 
states. 

The key feature of a thin superconducting film of thickness d is that the magnetic 
fields, as they live in 3 spatial dimensions, spread out in space outside the film and 
thus make the penetration depth enhanced into A e // = X 2 /d ^> A. For a thin film A e // 
can be really large (for d ~ 10 — 100 Angstroms, A e // ~ 0(lcm)), and the dynamics 
of vortices in the system is equivalent to that of a superfluid. It becomes possible to 
see a BKT phase transition in a superconductor. This enhancement was noted in [36] 
and for a nice review see [37]. Now for a single vortex that has a magnetic flux of a 
one flux quantum $ , the magnetic field is spread on a large area A ~ -^e//- Thus 
the magnitude of the magnetic field inside the vortex core is of the order B ~ ^o/^B 




which is practically zero. Thus, when modeling vortices in a thin superconducting 
film it is a good approximation to set the local magnetic field to zero. 

The specific system studied here is a relativistic superfluid, but the features of non- 
relativistic superfluid vortices are a useful guide for interpreting our results. Specifi- 
cally, we find that vortex configurations have distinct length scales characterizing the 
variations of the condensate, charge density, and superfluid density. Each of these 
length scales is characterized by a l/y/i//i c — 1 near the critical point. However, 
each of these scales generically has a different dependence on the chemical potential 
for larger /i as is seen in Fig. ([5]). 

As with holographic dark solitons [5116] the behavior of the charge density depletion 
fraction of the vortex core is very different for the < 0\ > and < O2 > condensates. 
For both dark solitons and vortices, the < 0\ > condensate has near 100% charge 
depletion in the core for low temperature (large /i). For < 2 > the core density 
depletion is much more modest, near 40%. This is consistent with the picture that 
the < 0\ > condensate is BEC-like, comprised of a point-like boson. Similarly, this 
is suggestive that the < O2 > condensate is more BCS-like, comprised of non-local 
Cooper pairs. 

In [T7J the impact that varying the type of superfluid across the BEC-BCS crossover 
has on vortices was explored. In that non-relativistic system, it was found that for 
BEC superfluids there is a single length scale characterizing the variations of the 
condensate. However, for a BCS superfluid the condensate profile was characterized 
by two length scales. In the vortex's tail the characteristic length scale is set by the 
size of the gap. Near the core, as the system locally approaches the normal phase, 
kp determines the characteristic scale. In Fig ([6]) we see some evidence that similar 
things happen for holographic superfluids. Specifically we see that there is really one 
length scale characterizing both the core and tail region for an < 0\ > type conden- 
sate. For < O2 > the core and tail are characterized by two distinct scales. However, 

9 A similar scaling of the magnetic field with the area was used in [8], but there the area A is kept 
finite and thus, the external magnetic field is nonzero. 
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the picture is not quite as clean as one might have hoped. In the non-relativistic 
system, the core length scale is smaller than for the tail. In the holographic < O2 > 
superfluid we find the opposite. If one accepts the basic picture that that the two 
types of holographic superfluids are of BEC and BCS types, on might view the core 
length scale as giving information about the Fermi surface (in the large N limit). 
Specifically, we would conclude that the Fermi momenta is vanishing in the large N 
limit. The vanishing of kf is consistent with the absence of Friedel oscillations in 
vortices and dark solitons [5j|6]. For other recent discussions of Fermi surfaces in 
holographic quantum liquids see [38H1U] . It would be interesting to compare their 
results to the behavior we see in the vortex cores. 

We also study the manner in which vortices allow one to estimate the critical 
velocity of the superfluid. Heuristically, this is due to the fact that the superfluid's 
local velocity increases beyond its critical velocity as one approaches the core. We can 
then estimate the critical velocity by the radius where the superfluid density starts 
to drop off. At low temperatures the Landau criterion sets an upper bound for the 
critical velocity of the superfluid in terms of the lowest sound mode. The fact that 
the 0\ superfluid's critical velocity comes closest to saturating the Landau criterion 
indicates that it is more likely to be a BEC-like superfluid (whose critical velocity is 
set by a sound mode). 

In addition to studying the vortices scales of variation, we confirmed that their 
energy cost over the homogeneous solution diverges logrithmically. While this is 
generically for any 2 + 1 dimensional vortex, the holographic reason is that the exte- 
rior of an AdS-Schwarzschild black hole effectively acts as a finite box in the radial 
direction. This is a nice confirmation that the bulk vortices give rise to superfluid 
vortices in the boundary theory as opposed to superconductor vortices. On general 
grounds one expects vortices in a superconductor to be finite energy excitations. 

In addition it appears that one can have a BKT transition at temperatures below 
the temperature of spontaneous symmetry breaking. To definitively resolve this issue 
on may need to include 1/N effects as well as gravitational backreaction. As in [5j[6], 
the absence of gravitational back reaction limits our ability to take the zero tem- 
perature limit of vortex configurations. Including gravitational backreaction would 
allow us to see many features unobstructed by thermal fluctuations. In addition, the 
backreacted geometry would present a very novel example of black hole hair. For 
recent work on black holes with inhomogeneous hair see [41J. 

Throughout this paper and in studies of dark solitons we have repeatedly seen 
features suggesting that the the choice of quantization in AdS^ corresponds to two 
types of superfluids. It would be worthwhile to probe this feature more closely. Ideally 
one would like to understand if there is any holographic analog of the "unitarity" limit 
in the BEC-BCS crossover. One way to probe this physics might be to vary the m? 
parameter in the bulk Lagrangian (studies in this direction were performed in [2"T j H2" j 
S3]). If we associate the operator scaling dimension with the type of superfluidity, it 
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is interesting to speculate that the unitarity limit would correspond to saturating the 
BF bound, m 2 = —9/4, in AdS±. One way that one might try to see this would be 
to study the way that the critical velocity varies with m 2 . In [T7] it was shown that 
the critical velocity was not a monotonic function across the BEC-BCS crossover. 
The peak value occurrs in the unitarity limit. Since the system is relativistic one 
might be able to go all the way from BCS to BEC and to relativistic Bose-Einstein 
condensation (RBEC) [H]. 

Finally, it would be very interesting to explore all of these questions in a setting 
more appropriate for laboratory experiments. Specifically we anticipate all of these 
ideas have analogs in non-relativistic gauge-gravity duals. It would be very interesting 
to study vortices in the backgrounds described in [4"5H4"8] 
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A Gauge fixing 

In this appendix we will outline gauge fixing necessary to obtain the equations in 
section [3J Because vortex solutions are cylindrically symmetric in the (p, 9) plane, 
we can take all gauge invariant quantities independent of the angle 9. Since we are 
interested in static solutions, we will also take all the fields to be independent of the 
time coordinated Choosing the gauge, A z = 0, and using dgR = (since R is gauge 

10 More general condition would have been to take all the gauge invariant quantities independent 
of t, but this leads to trouble discussed in [I] 
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invariant) equations ( IT9|) . (1201) and (121!) become 

= zd z {Ld z R) + + -<9 p (p<9 p (-)) 

2r 2:^ 2 p z 



-X-U 2 + /(^X) 2 + ^ ~f X)2 + (A p - 9 pX ) 2 ) (47) 
^ / P 



= -<9 p (pcU p ) + \d e d z A e - (-) 2 d zX (48) 
p p l z 

= /<9 Z 2 A + -d p (pd p A t ) + \d 2 A t - (-) 2 A t (49) 
p p z 

= d z {fd z A p ) + 4<9 2 A> - - (-) 2 (A> - ^x) (50) 

p l p z 

= d z (fd z A g ) + p<9 p (i«9 p A,) - p5 p (-<M p ) - (*) 2 (A - d eX ) (51) 

= d z (f(*) 2 d zX ) - - p dM^) 2 (A p - d pX )) 

-Ld e ((-)(A e -d eX )). (52) 
p z z 

Furthermore, to model a vortex we want the radial current to vanish, that is, A p — 
d pX = 0. 

We may use the residual gauge invariance and the requirements of regularity 
along the horizon and in the vortex core to simplify the remaining equations. Using 
do(dgx ~ Aq) = 0, equation ( 152"|) implies 

d z {{-ff{z)d zX ) = 0. (53) 

Which has the solution 

R 2 f{z) 

We can note that C is independent of 9 since = dg(d z x — A z ) = d$d zX . Regularity 
of x a t the horizon forces us to set C(p) = 0. 

Cylindrical symmetry implies that doF p g = and deFo z = 0. Substituting these 
in (1501) . gives 

d z (fd z A p ) = 0, (55) 

which has a solutions d z A p = ^jr^r ■ Again, regularity at the horizon requires D = 0. 
This means that A p = A p (6,p). Thus, by a gauge transformation we can set A p = 

n = d pX 

The symmetry also imposes dgF z g = ded z A e = 0, This has the general solution 
A e = A ( d 1 \6, p) + A d 2 \z, p). From de(dex + A e ) = and d p x = we can conclude that 
A^p is independent of p and can be thus absorbed into X with a gauge transformation. 
Thus A e = A e (z,p). 
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Similarly, cylindrical symmetry also requires = dg(dgx + Ag) = dgX, which 
implies that x = a + flQ for a and ft constants. In order for $ to be a single valued 
function /3 has to be an integer n e Z, which is the winding number of the vortex. 
We will also set a = as a part of the gauge choice. Note that the field x( z -> x ) is 
completely determined by the equations of motion, and thus, cannot be thought of 
as a non-normalizable mode0 The final simplification is due to the observation that 
At — dox — A t is gauge invariant, and therefore is independent of 9 

Taking advantage of the cylindrical symmetry, the equations may be written 

= zd z (^d z R) + + -d p {pd p {-)) 

z l z l z p z 



_(*)( ig + I^ziO!) (56) 

z f p 2 



= fd 2 z A t + -d p (pd p A t ) - {-fA t (57) 
P z 

= d z (fd z A g ) + pd p (-d p A g ) - (-) 2 (Ag - n) . (58) 

p z 

After defining R = R/z, these equations are equivalent to (I22H24I) . 



B Boundary conditions for the % field 

According to the AdS/CFT dictionary the expectation value of the order parameter 
field is given by 

(OJ-e^R^x). (59) 

Thus, x(0, x) should be a quantity determined by the state one is looking at. Therefore 
one should not be able to choose the value of x(0, boundary condition. We 

will show in this appendix that this is indeed the case for a general time independent 
solution of the equations of motion. The equations of motion for x i n terms of the R 
field are in the gauge A z = given by 

d z (fR 2 d zX ) + di{R 2 {d t x - A t )) = 0. (60) 

Since fl60|) is elliptic, boundary conditions should be provided at every spatial bound- 
ary. Since at the horizon ( l60i) becomes parabolic, it provides an effective boundary 
condition at the horizon (this is usually referred as a regularity condition). Next we 
will show that the other boundary condition at the AdS boundary is fixed by the 
equations of motion for the O2 condensate and by a "regularity" condition for the 0\ 
condensate. 

11 In Appendix [Bl we show that for a general solution of the field equations IT9l) . (|20l) and (|2ip . 
the phase field x is completely determined by the equations of motion and there are no freedom to 
choose boundary conditions for \ a t the AdS boundary. 



25 



We know the asymptotic expansion for the fields, R = + zRp) + ... and 
= + zA^ + .... By regularity the expansion for \ should come in positive 

integer powers x — X(o) + z X(i) + •••• Plugging these expansions into f )60|) gives the 

following equation up to second order in the powers of z 

2R {l) R {2)Xm + 2R 2 W X(2) + diOR(i)(diX(o) - 4°)) = 0- (61) 

6 X(3)-R(i) + 2(-R^ 2 ) + 2-R(i)-R( 3 ))X(i) + 8R( 1 )R( 2 )X(2) 

+ 2d l (R il) R i2) (d lX (o) - 4 0) )) + d^Rf^dan) - = 0. (62) 

For the O2 condensate we have R(\)(x) = and thus ( IBTj) is trivially solved, while 
f l6^|) simplifies into X(i) = 0. This means that for O2, the equation of motion forces 
the Neumann boundary condition 

d zX {z = 0) = 0, (63) 

and thus there is no more freedom to give boundary conditions to the \ field at the 
AdS boundary. 

For the case of the 0\ condensate we are not able to derive such a boundary 
condition from the equations of motion. Rather, we will use a physical argument, 
which sounds quite sensible. We require that the flux of the electric U(l) current into 
the z direction vanishes as we approach the AdS boundary, since literally the space 
ends there. This can be thought of as requiring the total charge in the system to be 
conserved^ So we require 

lim / d 2 x^^J z = 0. (64) 

The bulk current is given by the expression J z = ^y^ = gg 22 (\E'*<9 2 \I/ — ^/d z ^/* — 
2iA z \^\ 2 ) = fR 2 d z x- Now (JMD implies that 

lim [ d 2 x^$-R 2 d zX = 0. (65) 
*>->° J z=zo z 2 

This can be true only if d z \{z = 0) = 0. 



C Large p behavior 

In this appendix we outline an alternate approach to determine the large radius behav- 
ior of vortex solutions. The main observation is that for large p the partial differential 
equations may be reduced to a set of coupled ordinary differential equations, which 
may be solve using Mathematica's NDSolve routine. 

12 A similar argument was used in |49j for the energy momentum tensor and conservation of the 
total energy, in order to derive the BF bound. 



26 



In more detail, the gauge fixed equations of motion are: 

d z (fd z R) + -d p (pd p R) -fz + {Ae ~ -*±)r=0 (66) 

p V p f J 

fd 2 z A t + l -d p {pd p A t ) - R 2 A t = (67) 

d z (fd z A e ) + pd p (~d p A e ) - R 2 (A e - n) = 0. (68) 

If we assume that that all of the fields become spatially homogeneous for large p, we 
see that the R and A t reduce to the equations solved by [2]. We can then develop an 
expansion around this solution, 

R = R° + 5 J^ + ..., A t = A°+ 5 -^ + ..., and A e = A° + 5 -^ + ... (69) 

pi p2 pi 

One would like to solve for Ag, 6R, SA t , and 5A e in terms of the homogeneous solu- 
tions. 

The equations become 

d z (fd z A° e ) - (R°) 2 (A° e - n) 

+ p\ (d z (fd z 6A e ) - (R°) 2 6A e - 2R°A° e 5R) + ... = (70) 
(A ) 2 

d z (fd z 6R)-(z-^jL)6R 
96040 

- (A° e - n) 2 R° + ... = (71) 
(fd 2 z 5A t - R 2 5A t ^ - 2R°A° t 8R + ... = 0. (72) 
Equating the coefficients of powers of p gives 

d z (fd z A° e )-(R°) 2 (A° e -n) = (73) 
d z {fd z 5A e ) - (R°) 2 5A e - 2R°A° 9 5R = (74) 

d z (fd z 5R) -{z- { -jJ-)5R + —7^^ - « ~ = (75) 

fd 2 5A t - R 2 5A t - 2R°A° t 5R = 0. (76) 

In addition to these differential equations one must also impose boundary and 
regularity conditions. For the scalar field we should impose the boundary conditions 
appropriate for the type of condensate, d z 5R(0) — (< Ox > case) or 5R(0) = 
(< O2 > case). For 5A t we require that the chemical potential not be spatially 
varying, 5A t (0) = 0. Regularity at the horizon requires 5A t {l) = 0. To study 
vacuum properties we should set all external sources for the superflow to zero, Ajj(0) = 
6A e (0) = 0. 
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The most important thing to note is that the partial differential equations reduce 
to a set of coupled ODE's, which may be solved using any standard numerical dif- 
ferential equations solver. We have checked that the solutions to (175|) determined 
using Mathematica's NDSolve agree with the solutions described in Section 14.11 The 
difference between solutions found with the two methods is at the .1% level (the level 
to which the equations are solved by Gauss-Seidel method). Because the asymptotic 
analysis is wholly independent of the Gauss-Seidel approach, the agreement between 
the two provides an independent check of the results in Section 14. 11 



28 



References 



S. S. Gubser, "Breaking an Abelian gauge symmetry near a black hole horizon," 
Phys. Rev. D 78, 065034 (2008) |arXiv:0801.2977l [hep-th]]. 

S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, "Building a Holographic 
Superconductor," Phys. Rev. Lett. 101, 031601 (2008) jarXiv:0803.3"295 [hep- 
th]]. 

C. P. Herzog, P. K. Kovtun and D. T. Son, "Holographic model of superfluidity," 
larXiv:0809.4870l [hep-th]. 

P. Basu, A. Mukherjee and H. H. Shieh, "Supercurrent: Vector Hair for an AdS 
Black Hole," Phys. Rev. D 79, 126004 (2009) [arXiv:0809.4494l [hep-th]]. 

V. Keranen, E. Keski-Vakkuri, S. Nowling and K. P. Yogendran, "Dark Solitons 
in Holographic Superfluids," Phys. Rev. D 80, 121901 (2009) [arXiv: 0906 .52171 
[hep-th]]. 

V. Keranen, E. Keski-Vakkuri, S. Nowling and K. P. Yogendran, "Inhomoge- 
neous Structures in Holographic Superfluids: I. Dark Solitons," larXiv:0911.1 866 
[hep-th] . 

T. Albash and C. V. Johnson, "Vortex and Droplet Engineering in Holographic 
Superconductors," larXiv:09 06.1795 [hep-th]. 

M. Montull, A. Pomarol and P. J. Silva, "The Holographic Superconductor 
Vortex," larXiv:0906.2"396l [hep-th]. 

K. Maeda, M. Natsuume and T. Okamura, "Vortex lattice for a holographic 
superconductor," larXiv:0910.4475l [hep-th]. 

T. Albash and C. V. Johnson, "A Holographic Superconductor in an External 
Magnetic Field," JHEP 0809 (2008) 121 [ arXiv: 0804 .34661 [hep-th]]. 

S. A. Hartnoll, C. P. Herzog and G. T. Horowitz, "Holographic Superconduc- 
tors," JHEP 0812 (2008) 015 [ arXiv:0810.1563l [hep-th]]. 

T. Albash and C. V. Johnson, "Phases of Holographic Superconductors in an 
External Magnetic Field," arXiv:0906.051~9l [hep-th], 

H. b. Zeng, Z. y. Fan and Z. z. Ren, "Time Reversal Symmetry Breaking Holo- 
graphic Superconductor in Constant External Magnetic Field," Phys. Rev. D 
80 (2009) 066001 |arXiv:0906.2323l [hep-th]]. 



29 



[14] M. H. Dehghani, A. M. Ghezelbash and R. B. Mann, "Vortex holography," 
Nucl. Phys. B 625 (2002) 389 |arXiv:hep-th/0105134|. 



[15] M. H. Dehghani, A. M. Ghezelbash and R. B. Mann, "Abelian Higgs 
hair for AdS-Schwarzschild black hole," Phys. Rev. D 65 (2002) 044010 
larXiv:hep-th/010722"4] . 

[16] M. Antezza, F. Dalfovo, L. P. Pitaevskii, and S. Stringari, "Dark solitons in a 
superfTuid Fermi gas," Phys. Rev. A 76, 043610 (2007) |arXiv:0706.060T] . 

[17] R. Sensarma, M. Randeria, and T. L. Ho, "Vortices in Superfluid Fermi 
Gases through the BEC to BCS Crossover," Phys. Rev. Lett. 96, 090403 (2006) 
| arX iv: cond-mat /0510761] . 

[18] I. Amado, M. Kaminski and K. Landsteiner, "Hydrodynamics of Holographic 
Superconductors," JHEP 0905, 021 (2009) [ arXiv:0903.2209l [hep-th]]. 

[19] J. M. Kosterlitz, D. J. Thouless, "Ordering, metastability and phase transitions 
in two-dimensional systems," Journal of Physics C: Solid State Physics, Vol. 6 
pages 1181-1203 (1973). 

[20] V. L. Berezinskii, "Destruction of Long-range Order in One-dimensional and 
Two-dimensional Systems having a Continuous Symmetry Group I. Classical 
Systems," Soviet Journal of Experimental and Theoretical Physics, Vol. 32, 
p.493. 

[21] G. T. Horowitz and M. M. Roberts, "Holographic Superconductors with Various 
Condensates," Phys. Rev. D 78 (2008) 126008 jarXiv:0810. 10771 [hep-th]]. 

[22] C. P. Herzog, "Lectures on Holographic Superfluidity and Superconductivity," 
J. Phys. A 42 (2009) 343001 [ arXiv:0904.1975l [hep-th]]. 

[23] S. A. Hartnoll, "Lectures on holographic methods for condensed matter 
physics," Class. Quant. Grav. 26, 224002 (2009) |arXiv:0903.3"246l [hep-th]]. 

[24] S. S. Gubser, C. P. Herzog, S. S. Pufu and T. Tesileanu, "Superconductors from 
Superstrings," Phys. Rev. Lett. 103 (2009) 141601 |arXiv:0907.3510l [hep-th]]. 

[25] J. P. Gauntlett, J. Sonner and T. Wiseman, "Holographic superconductivity in 
M-Theory," Phys. Rev. Lett. 103, 151601 (2009) |arXiv:0907.3796l [hep-th]]. 

[26] S. S. Gubser, S. S. Pufu and F. D. Rocha, "Quantum critical superconductors 
in string theory and M-theory," Phys. Lett. B 683 (2010) 201 |arXiv:0908.00~TT1 
[hep-th]]. 



30 



[27] J. Gauntlett, J. Sonner and T. Wiseman, "Quantum Criticality and Holographic 
Superconductors in M-theory," larXiv :0912.0512 [hep-th]. 

[28] F. Denef and S. A. Hartnoll, "Landscape of superconducting membranes," Phys. 
Rev. D 79, 126008 (2009) jarXiv: 0901 .11601 [hep-th]]. 

[29] D. T. Son, "Low-energy quantum effective action for relativistic superfluids," 
arXiv:hep-ph/0204199"l 

[30] M. A. Valle, "Hydrodynamic fluctuations in relativistic superfluid," Phys. Rev. 
D 77, 2, (2008) [arxiv:hep-ph/0707.2665]. 

[31] F. Dalfovo, S. Giorgini, L. Pitaevskii, S. Stringari "Theory of Bose- 
Einstein condensation in trapped gases," Rev. Mod. Phys. 71, (1999) 
[arxi jcond-mat/980603 8] . 

[32] G. H. Derrick, "Comments on Nonlinear Wave Equations as Models for Ele- 
mentary Particles," J. Math. Phys. 5, 1252 (1964). 

[33] R. B. Diener, R. Sensarma, M. Randeria, "Quantum fluctuations in the su- 
perfluid state of the BCS-BEC crossover," Phys. Rev. A 77, 023626 (2008) 
|arXiv:0 709.2653 [cond-mat. other]]. 

[34] C. P. Herzog and A. Yarom, "Sound modes in holographic superfluids," Phys. 
Rev. D 80 (2009) 106002 [ arXiv:0906.4810l [hep-th]]. 

[35] T. Nishioka, S. Ryu and T. Takayanagi, "Holographic Superconductor/Insulator 
Transition at Zero Temperature," larXiv:0911.0 962 [hep-th]. 

[36] J. Pearl, Appl. Phys. Lett. 5, 65 (1964). 

[37] D. R. Nelson, "Defects in Superfluids, Superconductors and Membranes," 
[arXiv:9502114vl [cond-mat]]. 

[38] J. W. Chen, Y. J. Kao and W. Y. Wen, "Peak-Dip-Hump from Holographic 
Superconductivity," larXiv:0911.282T1 [hep-th], 

[39] T. Faulkner, G. T. Horowitz, J. McGreevy, M. M. Roberts and D. Vegh, "Pho- 
toemission 'experiments' on holographic superconductors," arXi v:0911.3 402 
[hep-th] . 

[40] S. S. Gubser, F. D. Rocha and P. Talavera, "Normalizable fermion modes in a 
holographic superconductor," larXiv:0911.3 632 [hep-th]. 

[41] S. Nakamura, H. Ooguri and C. S. Park, "Gravity Dual of Spatially Modulated 
Phase," larXiv:0911.0679l [hep-th]. 



31 



O. C. Umeh, "Scanning the Parameter Space of Holographic Superconductors," 
JHEP 0908 (2009) 062 [arXiv:0907.3136l [hep-th]]. 

Y. Kim, Y. Ko and S. J. Sin, "Density driven symmetry breaking and But- 
terfly effect in holographic superconductors," Phys. Rev. D 80 (2009) 126017 
|arXiv:0904.4567l [hep-th]]. 

Y. Nishida and H. Abuki, "BCS-BEC crossover in relativistic superfluid 
and its possible realization in QCD," Phys. Rev. D 72 (2005) 096004 
|arXiv:hep-ph/0504083] . 

D. T. Son, "Toward an AdS/cold atoms correspondence: a geometric re- 
alization of the Schroedinger symmetry," Phys. Rev. D 78, 046003 (2008) 
|arXiv:0804.3972l [hep-th]] . 

K. Balasubramanian and J. McGreevy, "Gravity duals for non- relativistic 
CFTs," Phys. Rev. Lett. 101, 061601 (2008) [arX iv:0804.4053l [hep-th]]. 

S. Kachru, X. Liu and M. Mulligan, "Gravity Duals of Lifshitz-like Fixed 
Points," Phys. Rev. D 78, 106005 (2008) |arXiv:0808. 17251 [hep-th]]. 

E. J. Brynjolfsson, U. H. Danielsson, L. Thorlacius and T. Zingg, "Holographic 
Superconductors with Lifshitz Scaling," larXiv :0908.2611 [hep-th]. 

O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, "Large 
N field theories, string theory and gravity," Phys. Rept. 323, 183 (2000) 
|arXiv:hep-th/9905111| . 



32 



